00 

o 

O 
(N 



(N 






IISER/HEP/07/08 

Lepton Flavours at the Early LHC Experiments as the 
Footprints of the Dark Matter Producing Mechanisms 

Nabanita Bhattacharyya^"-'^, Amitava Datta*^")^ and Sujoy Poddar*^"-'^ 

("'' Indian Institute of Science Education and Research, Kolkata, 
HC-VII, Sector III, Salt Lake City, Kolkata 700 106, India. 

Abstract 

The niSUGRA parameter space corresponding to light sleptons well within the 
reach of LHC and relatively light squarks and gluinos (mass < 1 TeV) has three regions 



1-^ ' consistent with the WMAP data on dark matter relic density and direct mass bounds 

Oh: 

from LEP 2. Each region can lead to distinct leptonic signatures from squark-gluino 



events during the early LHC experiments (integrated luminosity ~ 10 fb or even 
smaller). In the much studied stau-LSP coannihilation region with a vanishing common 



^ ' trilinear coupling {Aq) at the GUT scale a large fraction of the final states contain 
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t^^ . electrons and / or muons and e - fi - t universality holds to a good approximation. In 
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^D ' the not so well studied scenarios with non-vanishing Aq both LSP pair annihilation and 



stau-LSP coannihilation could contribute significantly to the dark matter relic density 
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Q I for even smaller squark-gluino masses. Our simulations indicate that the corresponding 

signatures are final states rich in r-leptons while final states with electrons and muons 

r> I are suppressed leading to a violation of lepton universality. These features may be 

observed to a lesser extent even in the modified parameter space (with non-zero ^o) 

where the coannihilation process dominates. We also show that the generic ?7i-leptons + 

n-jets+ ^ signatures without flavour tagging can also discriminate among the three 

scenarios. However, the signals become more informative if the r and b-jet tagging 

facilities at the LHC experiments are utilized. 
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1 Introduction 

Models with sup ersymnietry( SUSY) [1] are interesting for a variety of theoretical and 
phenomenological reasons. A specially attractive feature of the minimal supersymmetric 
standard model (MSSM) with R-parity conservation is the presence of the stable, weakly 
interacting lightest neutralino (x? ) [2] which is assumed to the lightest supersymmetric 
particle (LSP). This turns out to be a very good candidate for the observed dark matter 
(DM) in the universe [3, 4, 5]. 

Various SUSY models have been proposed and constrained by the data on DM relic 
density [6]- [14]. The recent revival of interest in this field is due to the very restrictive data 
from the Wilkinson Microwave Anisotropy Probe (WMAP) observation [15]. Combining the 
WMAP data with the results from the SDSS (Sloan Digital Sky Survey) one obtains the 
conservative 3 a limits 

0.09 < VLDKih^ < 0.13 (1) 

where QuMh"^ is is the DM relic density in units of the critical density, h = 0.71 ± 0.026 
is the Hubble constant in units of 100 Km s~^ Mpc~^. In this paper we shall assume that 
Qdm = ^^0. We should note here that the upper bound on Q^o in Eq.(l) must hold in any 
model with SUSY. In contrast the lower bound evaporates if the possibility of non-SUSY 
origin of DM is left open. 

In the thermally generated DM scenario the present value of ^^^oK^ can be computed 
by solving the Boltzmann equation for n^^o, the number density of the LSP in a Friedmann- 
Robert son- Walker universe. The most important particle physics input in this calculation is 
the thermally averaged quantity < UeffV >, where v is the relative velocity between two neu- 
tralinos annihilating each other and a^jf is the annihilation cross-section for all possible final 
states involving SM particles only. In addition to the negation of a LSP pair, coannihilation 
of the LSP [6]-[9],[10] with supersymmetric particles (sparticles) approximately degenerate 
with the LSP may also be important. The smaller the annihilation/coannihilation cross- 
section the larger becomes the LSP relic density. In addition to the parameters of the 
standard model the annihilation cross-section a^ff depends on the masses and the couplings 
of the sparticles and on the magnitudes of the bino (5), wino {W) and Higgsino {H^ , if 2) 
components of the LSP. Discovery of SUSY at the LHC followed by the measurement of 
the above parameters can, therefore, verify the hypothesis of supersymmetric DM as well as 
identify the underlying DM relic density producing mechanism [16]. This, however, is likely 
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to take some time. 

SUSY models with relatively small squark, gluino masses are of considerable contempo- 
rary interest since such strongly interacting sparticles with large production cross-sections 
are expected to show up in the early stages of the LHC experiments. In this case there are 
a few important relic density producing mechanisms like the LSP pair annihilation, LSP 
- lighter stau (fi) and LSP-lighter-stop (ti) coannihilation [13], [6]-[9], [10]. Each mecha- 
nism is active in one or more regions of the SUSY parameter space. It is worthwhile to 
check whether the sparticle spectra corresponding to each region yield distinct signatures 
at the early LHC experiments so that some idea of the DM producing mechanism, though 
qualitative, can be obtained. 

In this paper we shall restrict ourselves to the popular minimal supergravity 
(mSUGRA)model [17] with moderate values of the parameter tan (3 (to be defined in the 
next section). Our attention will be focussed on the regions of the parameter space consistent 
with the WMAP data and within the reach of the early runs. This region necessarily corre- 
sponds to light sleptons well within the reach of LHC but unlikely to be discovered directly 
by the early experiments (see section 3 for the details). Moreover the slepton signature alone 
carries very little information about the underlying DM relic density producing mechanism. 

In mSUGRA the above scenario also involves relatively light squarks and gluinos. The 
lighter chargino (xi ) and the second lightest neutralino (^2) are copiously present in the 
squark-gluino decay cascades. In the light slepton scenario these inos almost exclusively 
decay leptonically via two body decay modes. The lepton flavour content of the final states 
thus obtained from squark-gluino events at the early stages of the LHC run may reflect the 
underlying relic density producing mechanisms [18] (see section 2 for a brief review). Thus 
at least the footprints of these mechanisms may be viewed long before reconstruction of the 
sparticle masses and other relevant parameters establish the model rigorously. 

In the next section we shall review the important DM relic density producing mechanisms 
in different regions of the mSUGRA parameter space with sparticle spectrum as described 
above and qualitatively review the characteristics of the squark-gluino signatures from each 
region. In section 3 we shall go beyond [18] and present the results of our simulations 
revealing new aspects of the signals. This will justify the qualitative discussions of the 
previous section. The summary along with future outlooks will be the content of the last 
section. 



2 The Early LHC signatures: a qualitative discussion 

The simplest gravity mediated SUSY breaking model - the minimal supergravity 
(mSUGRA) [17] model- has only five free parameters. These are mi/2 ( the common gaugino 
mass), rriQ (the common scalar mass) and the common trilinear coupling parameter Aq, all 
given at the gauge coupling unification scale {Mg ~ 2 x 10^^ GeV), the ratio of Higgs vac- 
uum expectation values at the electroweak scale namely tan/5 and the sign of /i, the higgsino 
mixing parameter. The magnitude of fi is determined by the radiative electroweak symme- 
try breaking (REWSB) condition [19]. The low energy sparticle spectra and couplings at 
the electroweak scale are generated by renormalization group evolutions (RGE) of the soft 
breaking masses and coupling parameters [20]. In this paper we shall restrict ourselves to a 
moderate value of tan (3 namely 10 and positive fi. Representative values of the remaining 
parameters will be used to highlight different DM relic density producing mechanisms and 
the corresponding collider signals. Since the entire sparticle spectra and the couplings can 
be computed in terms of the five parameters only, the calculation of the LSP annihilation 
or coannihilation cross-sections and, consequently, the DM relic density are rather precise in 
this framework [6, 13, 21]. 

The WMAP allowed regions of the mSUGRA parameter space can be classified into 
several regions depending on the dominant LSP annihilation/coannihilation mechanisms. 
The details can be found in [3, 4, 5, 21]. In the following we list the mechanisms which will 
be relevant for the discussions in this paper focussing on relatively light sleptons, squarks 
and gluinos. 

One such region corresponds to small m-o but somewhat larger mi/2 (numerical examples 
will be given later). This choice leads to sleptons lighter than the charginos and the heavier 
neutralinos. Moreover, the mass difference between the lighter stau (fi) and the LSP turns 
out to be at most 30 GeV or so. Consequently ti -LSP coannihilation [6] often abbreviated 
as f-coannihilation occurs quite efficiently sufficient for producing the DM in the universe. 
The importance of this region, which we shall refer to as the conventional r- coannihilation 
zone, in the context of the WMAP data has recently been emphasized by several groups [21]. 

Most of the above analyses, however, were based on the choice Aq = without any 
compelling theoretical or empirical reason. One consequence of this ad hoc choice is that 
in order to satisfy the LEP bound on the lightest Higgs boson mass: rrih > 114.4 GeV [22] 
one requires relatively large mi/2- Typically for niQ ~ 100 GeV one requires mi/2 ~ 500 



GeV to satisfy the rrih bound and the WMAP constraints. As a result the squarks and 
gluinos become approximately degenerate, the latter being slightly heavier. The masses of 
these superpartners turn out to be (9(1 TeV) or larger. Nevertheless the total squark-gluino 
production cross section is sufficient to produce observable signatures at the LHC (See, e.g., 
[23]) throughout the f-coannihilation region. The lighter top squark mass is smaller than 
the other squarks due to the usual renormalization group effects driven by its large Yukawa 
coupling. However, due to the choice Aq = its mass is not further suppressed by mixing 
effects in the mass matrix. Therefore its mass is also of the order of one TeV. Such heavy ti 
does not participate in ti - LSP coannihilation. Moreover the production of ti pairs do not 
affect the total squark-gluino production significantly. The features discussed in the last two 
paragraphs lead to distinct collider signatures as we shall illustrate below with numerical 
examples. 

The other region of interest is the bulk annihilation region or the bulk region [4, 5, 13] 
where rriQ and mi/2 are such that many of the sparticles are significantly lighter than those 
in the conventional f-coannihilation region. The LSP turns out to be bino dominated and, 
consequently, couples favourably to right sleptons, which in fact are the lightest sfermions in 
this region of the parameter space. As a result an LSP pair efficiently annihilates into SM 
fermions via the exchange of light sfermions in the t-channel. This cross-section depends on 
the mass of the LSP (m~o), the masses of the exchanged sfermions and the LSP-sfermion 
couplings [4, 5, 13]. This region characterized by relatively light sparticles is especially 
interesting since SUSY signals with large events rates may be expected at the early LHC 
runs. 

Strong lower bounds on sparticle masses [24], particularly on the slepton masses from 
LEP 2 disfavor a part of the bulk annihilation zone. We, however, wish to emphasize that 
the direct bounds on the slepton masses alone can not eliminate the entire bulk region. 
Nevertheless a more severe restriction practically rules out the (mo — ^^11/2) plane containing 
the bulk region for Ao = 0. This arises from the LEP 2 bound on lightest Higgs boson mass 
(m/i) [22] since rrih and slepton masses are correlated in mSUGRA. Thus it has often been 
claimed in the recent literature that the mSUGRA parameter space with low values of both 
mo and 1711/2 and, consequently, the entire bulk region is strongly disfavoured [21]. 

It was emphasized in [18] that the above conclusions are artifacts of the ad hoc choice 
Aq = 0. On the other hand it is well known that for given rriQ and mi/2 moderate to 



large negative'* values of Aq lead to larger rrih [25]^. Hence in this case the bound on rrih 
can be satisfied even for relatively small mo and mi/2- This revives the region where LSP 
pair annihilation is the dominant DM producing mechanism. This can be seen, e.g., from 
Fig.l of [18] (see the blue(deep shaded)region; several other figures of the above reference 
corresponding to different choices of mSUGRA parameters also exhibit similar features). 
Moreover the low mo — m,i/2 regions of the mSUGRA parameter space are characterized by 
relatively light squarks and gluinos. This along with the inevitable presence of light sleptons 
in this region of the parameter space leads to distinct signals from squark-gluino events. 

There is also a region where f-coannihilation is still the most important mechanism for 
creating the observed DM in the universe( see the pink (the light shaded) regions of Fig. 1 
[18] and other similar figures). Remarkably, even this region corresponds to much smaller 
m,i/2 compared to what one would obtain for the Aq =0 case (i.e., for the conventional 
f-coannihilation region) . 

Furthermore large negative values of Ao leads to a relatively light top squark. In fact 
for a small but non-negligible region of the parameter space, the LSP - ti coannihilation [10] 
along with bulk annihilation may significantly contribute to the observed DM density with 
the ti well within the kinematic reach of the Tevatron (see the red region of Fig. 1 [18] and 
other figures). 

In [18] some features of the sparticle spectrum and signals at the Tevatron and the LHC 
corresponding to the WMAP allowed regions of the parameter space opened up by non-zero 
Aq were studied. The results were compared and contrasted with the expectations from 
the well publicized conventional fi-coanihilation scenario with Aq = hj introducing three 
benchmark scenarios A, B and C ( Table 1 of [18], which has also been reproduced here as 
Table 1 for a ready reference). In scenario A with a relatively large Aq both bulk annihilation 
and f-coannihilation are responsible for producing the observed DM relic density although 
the former dominates. On the other hand in scenario B with somewhat larger ?ni/2 the latter 
dominates but the contribution of the bulk annihilation is non- negligible. Finally in scenario 
C with Aq = the conventional f- coannihilation is the only significant DM producing 
mechanism. The scenario C corresponds to the smallest mi/2 which is consistent with the 



''We follow the standard sign convention of Ref. [26] for the signs of /i and Aq. 

^For Aq > 0, one requires toq and rriiM typically larger than the corresponding values for ^o=0- This 

does not lead to any novel collider signal. 



Higgs mass bound and the WMAP data. 

In the following the main features of the signal in the three scenarios [18] are summarized. 
The tables referred to in the rest of this section belongs to the original work unless stated 
otherwise explicitly. 

The sparticle spectra in the three scenarios can be found in Table 2 and the total squark- 
gluino cross-section in Table 6. These lowest order cross-sections have been computed by 
CalcHEP (version 2.3.7) [27]. 

The signals at the LHC are governed by the cascade decays of the above sparticles. 
In all three cases the gluinos are heavier than all squarks. As a result the gluinos decay 
into quark-squark pairs (see Table 3) . The squarks belonging to the third generation are 
relatively light due to the usual renormalization group effects in the mSUGRA model. Their 
masses are further suppressed by mixing effects in the mass matrix due to the non-vanishing 
Aq parameter scenarios A and B. As a result these squarks are more frequently present in 
gluino decay products in these scenarios compared to C 

The squarks belonging to the first two generations in general decay into the corresponding 
lighter quarks and an appropriate electroweak gaugino (Table 3). The lighter top squark 
similarly decays into appropriate quark-gaugino pairs. The bottom squarks of both type 
may decay, in addition to above channels, into a lighter top squark and a W boson in some 
scenarios (Table 4) with non- negligible branching ratios(BRs). Nevertheless, the decay of 
each third generation squark inevitably contains a b quark. This is the origin of the large 
fraction of final states with b-partons as noted in [18]. In scenario C the fraction of third 
generation squarks in gluino decay is relatively small and the above effect is suppressed to 
some extent. 

The decay properties of the lighter chargino (xi^)and the second lightest neutralino (x^) 
which, in addition to the LSP, are often present in squark-gluino decays determine the lepton 
content of the final sates to a large extent. As a direct consequence of the presence of the 
light sleptons, these two unstable gauginos decay almost exclusively into leptonic channels 
via two body modes in all three scenarios. 

In Table 2 of this paper, which is an enlarged version of Table 5 of [18], we present the 
relevant branching ratios of Xi a-nd X2- ^^ scenario A the wino dominated lighter chargino, 
is not kinematically allowed to decay into L-type charged sleptons. Its two body decay into a 
lepton-sneutrino pair, though kinematically allowed, is phase space suppressed. It therefore 



decays into R-type sleptons through its sub dominant Higgsino component with a large BR. 
This resuhs in a very large fraction of final states containing the ri, which eventually decays 
into a r-LSP pair. Only a tiny fraction of the final states contains electrons and muons. 
Unlike the SM the hallmark of this scenario is, therefore, lepton non-universality in the 
final states. For similar reasons the X2 decays primarily into t-ti pair while a much smaller 
fraction decays into neutrino-sneutrino pairs, contributing further to the r dominance of the 
final states. The sneutrinos in turn decay into the invisible neutrino-LSP channel in all three 
scenarios and act as additional carriers of missing energy. 

The scenario B has all the above features albeit to a lesser extent. The r dominance 
in both Xi ^-nd X2 decays exists but is reduced significantly compared to the predictions of 
scenario A. The fraction of invisible decays of X2 i'^^o v — v pair increases. This common 
feature of A and B can be easily illustrated by a parton level calculation (see Table 7). In a 
realistic LHC experiment with good r tagging capabilities the observability of this r excess 
has been demonstrated by simulation(see Table 8) using Pythia. 

In Table 2 our focus was on the three bench mark scenarios. In order to illustrate that 
the lepton flavour content of the final state is indeed correlated with the DM producing 
mechanism we have selected several representative points [S\ — S^) from the figures in [18]. 
The dominant DM producing mechanism and the figure no. in [18] are given in parentheses. 
For all points tan/5 = 10; /i > 0. The other parameters are as given below. 
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mo = 100; mi/2 = 250; Aq = -700 (Bulk ; Fig. 4) 

rriQ = 125; mi/2 = 400; Aq = -700 (Coannihilation ; Fig. 4) 

rriQ = 120;mi/2 = 400; Aq = -800 (Coannihilation ; Fig. 1(a)) 

mo = 120;mi/2 = 265; Aq = -900 (Bulk ; Fig. 1(a)) 

mo = 120;mi/2 = 300; Ao = -1000 (Bulk ; Fig. 1(a)) 

mo = 170;mi/2 = 500; Ao = -1200 (Coannihilation ; Fig. 1(c)) 

mo = 170;mi/2 = 500; Ao = -1400 (Coannihilation ; Fig. 1(c)) 



The relevant BRs of lighter chargino and second lightest neutralino decays are presented 
in Table 3. It is readily seen that if the bulk annihilation (f -coannihilation ) is the dominant 
mechanism, the decays follow the patterns of scenario A (B) (compare with Table 2). 

In scenario C the xf decays into left slepton- neutrino pairs or sneutrino -lepton pairs 



belonging to the first two generations with equal BR of sizable magnitudes. The fraction 
of final states involving rs is only marginally larger and lepton universality holds to a very 
good approximation. In fact due to limited r detection efficiency, the fraction of observed 
final states involving e and/or /i will be apparently larger in stark contrast to the predictions 
of scenarios A and B. The decays of X2 contributes further to the restoration of lepton 
universality. The X2 ^ow dominantly decays into the invisible final state consisting of a 
neutrino and a sneutrino with 51.6 % BR. Thus it largely acts as a carrier of missing energy 
in addition to the LSP and the sneutrino. Although we have presented numerical results 
for scenario C only we have verified that the above features hold qualitatively for the entire 
r-coannihilation strip. 

Scenarios with the electroweak gauginos decaying into purely leptonic two body channels 
and sneutrinos and X2 acting as additional carriers of missing energy have already been 
discussed in details in the context of MSSM. The characteristic signatures from final states 
containing excess of electrons and /or muons and more than one invisible particles at the 
Tevatron [28] and at e"*" — e~ colliders like LEP or the NLC [29, 30] were simulated at the 
parton level. An ISAJET based analysis in the context of the Tevatron Run I was also done 
[31]. The scope of accommodating such scenarios in the mSUGRA models and several of its 
variants was also discussed [32]. 

However, the connection of this scenario with a WMAP allowed region of the mSUGRA 
parameter space has not been highlighted in the existing literature. Nor were the LHC 
signatures of this scenario studied with due emphasis. It is interesting to note that this virtual 
LSP('VLSP') or effective LSP('ELSP') scenario^ is realized in the well studied conventional 
f-coannihilation region of the mSUGRA parameter space. In the next section we shall 
simulate the LHC signatures of this scenario in great details. 

In this paper we complement and extend the analysis in [18] in several ways. First of 
all we have simulated the popular ?n-leptons + ri-jets + ^t signature for n > 2 and several 
choices of m, which stands for the number of electrons and muons. No fiavour tagging 
is required at this stage. It is gratifying to note that the scenarios can be distinguished 
reasonably well via these signatures inspite of the uncertainties in the cross-sections due to 
the choice of QCD scales. The relative enhancement of the final states with e and/ or fi 



^The invisibly decaying i/'s and X2 acting as additional carriers of missing energy may be called VLSP or 
ELSP in the context of collider physics 



compared to m = final states, a distinct characteristic of tlie VLSP scenario, lias been 
illustrated in the next section with several examples in scenario C 

Next we consider events of the type It + X, Ifi + X and le + X, where X includes all 
possible states with two or more jets but no stable lepton or tagged r ( the difference in the 
definition of X compared to that in [18] should be noted). In this work we also consider 
final states with 2 leptons of the same flavour in all possible charge combinations + X. We 
have used the r-tagging efficiencies provided by the CMS collaboration [33] in our Pythia 
based simulations. As discussed in the earlier paragraphs an excess of events with r have 
been demonstrated in scenario A and to a lesser extent in scenario B. Finally we examine 
the 6-jet content of the final states beyond the parton level using Pythia. 

In [18] the ti events were assumed to be the dominant source of backgrounds. In this 
paper we have extended the background analysis by considering several other processes. In 
particular we have found that for some signals the QCD background is more important than 
the ti background. The details will be presented in the next section. 

3 The Signals at the LHC 



mSUGRA 
parameters 


A 


B 


C 


mo 


120.0 


120.0 


120.0 


mi/2 


300.0 


350.0 


500.0 


Ao 


-930.0 


-930.0 


0.0 


tan/3 


10.0 


10.0 


10.0 


signifi) 


1.0 


1.0 


1.0 



Table 1: Three bench mark scenarios introduced in [18]. 



In this section we begin with the generic SUSY signals of the type m — I + n — j+ ^t 
,where / = e or /i and j is any jet. At first we do not employ any flavour tagging. Our aim is 
to study the feasibility of discriminating among the three models under consideration using 
these generic signals. Next we shall employ flavour tagging and demonstrate that it further 



Decay modes 
(Gauginos) 


A 


B 


C 


xl - i-J^ 


0.0 


2.2 


23.1 


xl - I-rI^ 


0.9 


0.6 


0.8 


xl - r^r- 


78.6 


46.0 


8.4 


Xl - r+r" 


0.0 


0.7 


10.2 


xl - x'lZ 


0.3 


0.2 


0.4 


xl - x\h 


0.0 


1.9 


5.2 


xl -^ ^m 


5.2 


24.0 


33.8 


xl -^ ^r^^r 


15.0 


24.0 


17.8 


xt - x\w+ 


2.6 


2.6 


5.1 


xt ^ ^//+ 


5.4 


25.0 


36.0 


xt -^ ^tT^ 


16.0 


25.0 


19.2 


xt ^ It^i 


0.0 


2.0 


22.0 


xt -^ T'l^r 


76.0 


44.0 


7.8 


xt -^ ^2^T 


0.0 


0.7 


9.6 



Table 2: The BRs of the dominant decay modes of the hghter chargino and the second 
hghtest neutrahno. All sneutrinos decay into the invisible channel v + ^in the three cases 
understudy. Here / stands both for e and /x. 

enhances our discriminatory power. 

The production and decay of all squark-gluino pairs are generated by Pythia (version 
6.409) [34]. Initial and final state radiation, decay, hadronization, fragmentation and jet 
formation are implemented following the standard procedures in Pythia. We have used the 
toy calorimeter simulation (PYCELL) provided in Pythia with the following criteria: 

• The calorimeter coverage is |?7| < 4.5. The segmentation is given by Arj x A0 = 
0.09 X 0.09 which resembles a generic LHC detector. 



A cone algorithm with A R= v^Arp~+"A0^ = 0.5 has been used for jet finding. 
ET*min = 30GeV and jets are ordered in Ex. 
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Decay modes 
(Gauginos) 


Si 


^2 


53 


54 


55 


56 


57 


xl - \i^ 


0.0 


10.2 


12.5 


0.0 


0.0 


15.5 


15.0 


xl - I'rI^ 


1.6 


0.7 


0.5 


0.9 


0.7 


0.2 


0.2 


xl - rtr- 


85.8 


26.6 


24.2 


95.4 


76.9 


19.7 


20.2 


Xl - r+r- 


0.0 


4.3 


5.4 


0.0 


0.0 


9.8 


10.6 


xl - x\z 


0.0 


0.2 


0.2 


0.2 


0.2 


0.1 


0.1 


xl - X'lh 


0.0 


3.3 


2.3 


0.0 


0.0 


2.0 


1.4 


xl -^ ^m 


2.8 


31.8 


31.3 


0.0 


5.5 


28.8 


27.3 


xl -^ r^rl^r 


9.7 


22.2 


22.7 


3.5 


16.7 


23.8 


25.0 


xt ^ xlw+ 


1.3 


3.6 


2.5 


1.9 


2.2 


2.0 


1.5 


xt ^ i^il^ 


3.0 


33.4 


33.4 


0.0 


5.7 


29.6 


28.0 


xt - ^.r+ 


11.0 


23.4 


23.9 


3.6 


17.7 


24.5 


25.7 


xt ^ /l^/ 


0.0 


10.1 


11.8 


0.0 


0.0 


15.1 


14.6 


X^ -^ rt^r 


84.7 


25.4 


23.1 


94.5 


74.3 


19.2 


19.8 


xt -^ ^2+//^ 


0.0 


4.0 


5.2 


0.0 


0.0 


9.5 


10.3 



Table 3: Same as Table 2 in the scenarios 5i — 57(see text). 

The stable leptons are selected according to the criterion : 

• Leptons (/ = e, /i) are selected with Pt > 30 GeV and |?7| < 2.5. For lepton-jet isolation 
we require AR{l,j) > 0.5. The detection efficiency of the leptons are assumed to be 
100%. 

The following cuts are implemented for background rejection : 

• Leptons (/ = e,yu) with Pt < 60 GeV are rejected to ensure the rejection of leptons 
coming from r decay. (CUT 1) 

• We reject events without at least two jets having Pt > 150 GeV( CUT 2) 

• Events with missing energy (^t) < 200 GeV are rejected. (CUT 3) 
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Events with Mejf < 1000 are rejected, where Me// = \^\ + Tji\P^\ + Si|P^'| (/ = e,/i 
).(CUT 4) 

Only events with jets having St > 2.0, where St is a standard function of the eigen- 
values of the transverse sphericity tensor, are accepted. (CUT 5) 



QCD scale 


0.5v^ 


V~s 


2.075 


A 


19.41 


15.58 


12.10 


B 


6.79 


5.74 


4.26 


C 


0.84 


0.74 


0.51 



Table 4: Variation of total cross-section in pb of squrak-gluino events with the QCD scale 
in A, B and C. 



Q 


0.5V§ 


v^ 


2.0v^ 


A 


0/ 


2.54 


2.03 


1.59 




1/ 


0.25 


0.20 


0.15 


B 


0/ 


1.45 


1.21 


0.90 




11 


0.21 


0.17 


0.13 


C 


01 


0.24 


0.21 


0.15 




11 


0.11 


0.09 


0.06 



Table 5: The cross-sections (including efficiency) in pb of events with 771 = and 1 (see 
text). 

Cut 1 will be employed in the second part of our analysis when the relative abundance 
of final states with e, fi and r will be studied. For establishing m — I + n — j+ ^t signals 
Cut 2 - Cut 4 are adequate. 

Table 4 illustrates that, as expected, the variation of the total cross-section of squark- 
gluino events with the choice of the QCD scale is rather large in each model. In spite of this 
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large variation it is clear that each scenario is characterized by a typical size of the cross- 
section. This cross-section is the largest in scenario A because squark and gluino masses are 
the smallest in this case (see [18]). In particular the contribution of the lighter top squark 
enhances the cross-section significantly. The last feature is a direct consequence of large 
negative Aq in scenario A. The corresponding masses are significantly larger in the other two 
scenarios resulting in smaller cross-sections. 

It is, however, impossible to conclusively identify a particular DM scenario by the size of 
the cross-section alone since a similar cross-section may arise from a different combination 
of mSUGRA parameters which may or may not be allowed by the relic density constraint. 
Moreover, the signal corresponding to a larger raw cross-section may eventually be suppressed 
due to the effects of the kinematical cuts, small BRs of the underlying decays etc. Several 
examples of this will be presented in the following paragraphs. 

The total cross-section at best provides a hint for the underlying SUSY model and LSP 
annihilation mechanism but no definite conclusion can be drawn. In the following we shall 
show that a multi-channel analysis using signals with different choices of m (the number of 
leptons in the final state) may very efficiently discriminate among different scenarios. 

First we assume that the SM backgrounds will be determined either from data or from 
theory with next to leading order accuracy and can be subtracted from the event sample 
without introducing large errors. The leading order signal cross-sections in pb (including 
the BRs and the efficiency of the cuts C2 -C5 listed above) for m = and 1 are presented 
in Table 5. It is readily seen that the variation of the cross-sections due to the QCD scale, 
the difference in BRs (see section 2 and [17] for further details) and the efficiency of the 
kinematical cuts in different scenarios may conspire in such a way that the cross-sections 
for m = (ctq) may look rather similar in different scenarios (see, e.g., o"o in model A at 
Q = 2.0vs and in model B at Q = 0.5Vs)at least for relatively low integrated luminosities 
(£). Similarly the predictions for m = l(cri) may appear to be quite similar in different 
scenarios. However, it can be easily checked that in each scenario the ratio R=(To/ci is scale 
independent to a very good approximation. The value of R for C of 10 fb~^ in A, B and 
C are 10.3, 7.0 and 2.3 respectively. Introducing a statistical uncertainty ^/N for N counts 
and using the standard method for estimating the uncertainty we find R^ = 10.3 ± 0.24, 
Kb = 7.0 ± 0.18 and Re = 2.3 ± 0.9. Thus it is fair to conclude that different scenarios can 
be readily distinguished from each other. Ratios involving cross-sections with other choices 
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of m exhibit similar scale independence. 

The relatively large value of R in scenario A is partly due to the r dominance of the final 
states in this scenario as discussed in the last section. Since the r decays into hadrons with 
a large BR, the numerator of R is naturally enhanced. The denominator of R on the other 
hand is small because, as explained in the last section, the number of final states involving 
e and /i are suppressed. 

The above properties also hold in scenario B to a lesser extent yielding a value of R 
smaller than that in A but significantly larger than the one in C In C, as in any other VLSP 
scenario, the denominator of R is rather large for reasons already discussed and a smaller R 
is obtained. 

In order to provide some estimate of the dominant backgrounds we present in Table 6 the 
signal and the important standard model backgrounds for several values of m in the leading 
order for Q = Vs. This will be followed by the usual analysis of the significance (S/yB) , 
where S{B) is the total number of signal (background) events. Although we have listed only 
the backgrounds from tt, QCD (including all quark-anti-quark and gluon events in Pythia 
in the lowest order ) we have also simulated WW, ZZ, WZ and Drell-Yan backgrounds and 
have found them to be indeed negligible. For later use we have simulated the background 
from W + jets events. 

The W + jets cross-section has been computed for pT > 50 where pT is defined in the 
rest frame of the parton -parton collision. The QCD cross-section has been computed in 
two Pt bins : (i) 400 GeV < pr < 1000 GeV and (ii) 1000 GeV< pr < 2000 GeV .The 
corresponding cross-sections being 2090 pb and 10 pb respectively. Outside these bins the 
number of events are negligible. 

Although the squark-gluino production cross-section is rather tiny in scenario C the 
signal cross-sections predicted for ?72 > 2 is larger than the corresponding signals in A and B 
with much larger raw production cross-section. This again is a direct consequence of the large 
leptonic BRs of gaugino decays in the underlying VLSP scenario as discussed above. The 
significance (without systematic errors) of the signal for different m and the representative 
£ of 10 fb^^ are in Table 7. The corresponding numbers for other C s can be estimated 
by simple scaling and one can easily verify that in several channels statistically significant 
signals can be obtained for much lower C s. 

From Tables 6 and 7 it is clear that one way to unambiguously discriminate between A, 
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B on the one hand and C on the other is the count of — Z events. This conclusion based 
on leading order cross-sections is likely to hold inspite of the scale uncertainty discussed 
above and the possibility that the systematic errors, which we have not considered in this 
paper, might be relatively large in the early stages of the LHC and affect the significance. 
The same count may discriminate between A and B although the theoretical uncertainties 
illustrated in Table 4 may cast some doubt on the results. The observation of the clean almost 
background free signal for m = 3 may vindicate model C since no statistically significant 
signal is expected from A or B in this case. On the other hand in the VLSP scenario C an 
acceptable signal in this channel is expected even for C significantly smaller than 10 fb~^. 

Counting experiments alone for m = 1 ot m = 2 may not be very useful for discrimi- 
nating between A and B during the early stages of the LHC run because of the theoretical 
uncertainties and low statistics. As discussed above one can form several QCD scale invariant 
ratios of observables to distinguish between scenarios A and B. For example R' = (Jq/(J2os 
,where (T20S is the cross-sections for m = 2 involving opposite sign leptons, R^ = 312.3 ±38.7 
and R'^ = 131.5 ± 13.8. Another example is R" = ai/a2ss ,where a2ss is the cross-sections 
for m = 2 same sign leptons, R^ = 83.3 ±17.1 and R'^ = 50.0 ±8.7. However, more statistics 
will be required to make the distinction unambiguous. 

One could use the next to leading order cross-section for squark-gluino production [35]. 
However, the uncertainty in the dominant QCD background in the canonical — / channel 
due to higher order effects is not known precisely. We have therefore restricted ourselves 
to leading order cross-sections for both the signal and the backgrounds. It is somewhat 
reassuring to note that the leading order signal cross-sections are typically multiplied by a 
factor of 1.4 - 1.5 in the next to leading order. Thus the significance will remain the same 
even if we have underestimated the background by a factor of two. 

Obviously the dominance of final states involving r leptons in some of the models under 
consideration can not be directly established from the generic observables. Nor can the 
dominance of final states with B-hadrons in certain scenarios be tested. For this r and 6-jet 
tagging facilities must be relied upon. 

We, therefore, turn our attention to final states of the type lr±X where X includes two 
or more hard jets but no e or /x or tagged r. Tagging of r jets are implemented according to 
the following procedure. 

Only hadronic r decays are selected. The r-jets with 77 <3.0 are then divided into 
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SIGNAL 


background 




A 


B 


C 


ti 


QCD 


a(pb) 


15.58 


5.74 


0.74 


400 


2100 


0/ 


2.03 


1.21 


0.21 


0.33 


3.55 


1/ 


0.20 


0.17 


0.09 


0.16 


5.0 X 10-3 


SS 


2.4 X 10-3 


3.4 X 10-3 


6.0 X 10-3 


7.2 X 10-^ 


- 


OS 


6.5 X 10"=^ 


9.2 X 10-3 


0.02 


0.01 


- 


2,1 


2.6 X 10"^ 


4.9 X 10-^ 


3.5 X 10-3 


3.2 X 10-^ 


- 



Table 6: The cross-sections (including efficiency) at Q = Vs for signal process with different 
m, ti and QCD events. Here SS refers to ?7i = 2 with leptons carrying the same charge and 
05* refers to similar events with leptons carrying opposite charge. No entry in a particular 
column (-) means negligible background. 

several Pt bins. A r-jet in any Pt bin is then treated as tagged or untagged according to the 
efficiency [e-r] given in [33] figure 12.9 for a particular bin. In this analysis we have applied 
Cuts 1-5. The corresponding efficiencies for various signals in scenario A and the dominant 
QCD and ti backgrounds are listed in Tables 8-10. 

The computation of le + X type events are rather straight forward. Here for simplicity 
we have assumed the e-detection efficiency to be 100 %. In our generator level analysis the 
result for 1/i + X is expected to be the same to a good approximation and we do not present 
them separately. The number of events of the above two types subject to the kinematical 
cuts listed above are presented in the second and third rows of Table 11 along with the 
dominant SM backgrounds. The QCD background to It + X events stems from mistagging 
of light flavour jets as r-jets. The mistagging probability has also been taken from [33] figure 
12.9. 

The It + X signal, if unambiguously observed, will disfavour model C. The size of 
the signal or the ratio N(lr + X)/N(le + X) can distinguish between scenarios A and B 
inspite of theoretical and possible systematic uncertainties. If this signal is not observed 
then the le + X signal can establish scenario C. If we require X to be free of 6-jets then 
the significance marginally increases because of the reduced background from tt events. The 



16 





A 


B 


C 


0/ 


103 


61 


11 


1/ 


49 


42 


22 


ss 


9 


13 


22 


OS 


5 


8 


16 


31 


1 


3 


20 



Table 7: The significance {^/yB{tt + QCD))oi signals in Table 6 for C= lOfb ^ . 





Efficiencies for selection cuts 


Cut 1 


Cut 2 


Cut 3 


Cut 4 


Cut 5 


Ir + X 


0.0712 


0.0441 


0.0334 


0.0263 


0.0155 


le + X 


0.0214 


0.0084 


0.0061 


0.0051 


0.0034 


lr + 06 + X 


0.0391 


0.0273 


0.0201 


0.0161 


0.0097 


iT+lh + X 


0.0194 


0.0112 


0.0081 


0.0062 


0.0034 


lT + 2h + X 


0.0113 


0.0061 


0.0041 


0.0039 


0.0022 


le + 06 + X 


0.0065 


0.0030 


0.0023 


0.0019 


0.0012 


le+16 + X 


0.0083 


0.0031 


0.0022 


0.0019 


0.0012 


le + 26 + X 


0.0047 


0.0019 


0.0014 


0.0012 


0.0008 



Table 8: The cumulative efficiency of the cuts for signal(A) given by Ni/N, where Xj is the 
number of events survived after successive application of Cut 1 to Cut i and X is the total 
sample generated. 



17 





Efficiencies for selection cuts 


Cut 1 


Cut 2 


Cut 3 


Cut 4 


Cut 5 


Ir + X 


0.1374 


0.0135 


0.0006 


0.0006 


1.29x10-^ 


le + X 


3.7x10-5 


3.6x10-5 


2.7x10-6 


2.7x10-6 


1.9x10-7 


lr + 06 + X 


0.1261 


0.0125 


2.7x10-^ 


2.4x10-^ 


4.6x10-5 


iT+lh + X 


0.0087 


0.0085 


2.8x10-^ 


2.5x10-^ 


5.6x10-5 


lr + 2h + X 


0.0019 


0.0019 


8.5x10-5 


7.9x10-5 


2.5x10-5 


le + 06 + X 


2.2x10^5 


2.1x10-5 


1.3x10-6 


1.3x10-6 


1.9x10-7 


le + lh + X 


1.2x10-5 


1.2x10-5 


1.3x10-6 


1.3x10-6 


- 


le + 2b + X 


2.0x10-6 


2.0x10-6 


- 


- 


- 



Table 9: Same as Table 8 for QCD events. 





Efficiencies for selection cuts 


Cut 1 


Cut 2 


Cut 3 


Cut 4 


Cut 5 


iT + X 


0.0392 


0.0041 


0.0006 


0.0005 


0.0001 


le + X 


0.0481 


0.0034 


0.0004 


0.0003 


0.0001 


iT + Ob + X 


0.0101 


0.0008 


0.0002 


0.0001 


2.3x10-5 


Ir + lb + X 


0.0192 


0.0023 


0.0003 


0.0002 


5.6x10-5 


lT + 2b + X 


0.0097 


0.0013 


0.0002 


0.0001 


3.0x10-5 


le + Ob + X 


0.0133 


0.0006 


0.0001 


6x10-5 


2.5x10-5 


le+16 + X 


0.0244 


0.0016 


0.0002 


0.0001 


4.6x10-5 


le + 26 + X 


0.0113 


0.0011 


0.0001 


9x10-5 


3.0x10-5 



Table 10: Same as Table 8 for tt events. 
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SIGNAL 


BACKGROUNDS 


A 


B 


C 


tt 


W + jets 


QCD 


a(pb) 


15.58 


5.74 


0.74 


400 


650 


2100 


It + X 


0.2415 


0.1251 


0.0232 


0.0433 


0.0018 


0.2709 


le + X 


0.0530 


0.0433 


0.0302 


0.0403 


0.0015 


0.0004 


Ir + 06 + X 


0.1511 


0.0744 


0.0151 


0.0091 


0.0014 


0.0958 


iT + lh + X 


0.0530 


0.0293 


0.0041 


0.0222 


0.0003 


0.1184 


lT + 2h + X 


0.0343 


0.0212 


0.0032 


0.0120 


0.0001 


0.0524 


le + 06 + X 


0.0187 


0.0141 


0.0201 


0.0101 


0.0013 


0.0004 


le+16 + X 


0.0187 


0.0152 


0.0052 


0.0182 


0.0002 


- 


le + 26 + X 


0.0125 


0.0121 


0.0043 


0.0120 


- 


- 



Table 11: The cross- sect ions (including efficiency) of events with one detected r and one 
isolated e. Here X stands for all possible final states excluding any lepton or tagged r but 
with at least two jets. The number of tagged 6-jets is given by n-6, n = 1, 2, 3 (see the ffist 
column of row 4 - 9) . 

statistical significance of various signals for the representative C of 10 fb^^ are listed in Table 
12. 

We next illustrate that both scenarios A and B are expected to be rich in 6-jets. A jet 
with |?7| < 2.5 matching with a B-hadron of decay length > 0.9mm has been marked tagged. 
The above criteria ensures that eh — 0.5 in ti events, where ej, is the single 6-jet tagging 
efficiency (i.e., the ratio of the number of tagged 6-jets and the number of taggable 6-jets in 
tt events). It is readily seen that, the fraction of le + X events with at least one tagged 6-jet 
is quite large in A (0.59) and B (0.63) and somewhat smaller in C(0.31). The significance of 
different signals are in Table 12. 

Finally we present in Table 13 events of the type 2r + X and 2e + X with different 
number of tagged 6-jets. It is to be noted that 2r + X events (2e + X events) are statistically 
significant in A and B (C). Moreover an observable signal in 2r +(> 16) +X channel is 
expected only in model A. 

Although all the scenarios under consideration involve light sleptons direct slepton 
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SIGNAL 


A 


B 


C 


a(pb) 


15.58 


5.74 


0.74 


Ir + X 


42.9 


22.2 


4.1 


le + X 


26.3 


21.0 


14.6 


lr + 06 + X 


46.4 


22.7 


4.6 


lr+(> lb)+X 


19.2 


11.0 


1.5 


le + Ob + X 


18.3 


12.8 


18.3 


le+(> 16) +X 


17.9 


15.6 


5.4 



Table 12: The S/VB{tt + Wjets + QCD) ratio for the signals in Table 11 corresponding to 





SIGNAL 


BACKGROUNDS 


A 


B 


c 


ti 


QCD 


a(pb) 


15.58 


5.74 


0.74 


400 


2100 


2t + X 


0.0171 


0.0078 


0.0019 


0.0015 


0.0231 


2e + X 


0.0008 


0.0009 


0.0031 


0.0010 


- 


2T + 0b + X 


0.0105 


0.0039 


0.0012 


0.0002 


0.0084 


2T + lb + X 


0.0038 


0.0019 


0.0003 


0.0007 


0.0073 


2T + 2b + X 


0.0028 


0.0017 


0.0003 


0.0006 


0.0067 


2e + 0b + X 


0.0004 


0.0004 


0.0021 


0.0002 


- 


2e+lb + X 


0.0003 


0.0003 


0.0005 


0.0006 


- 


2e + 2b + X 


0.0001 


0.0002 


0.0004 


0.0002 


- 



Table 13: The cross-sections (including efficiency)of events with two detected r and two 
isolated e. The other conventions are as in Table 11. 

searches are unlikely to yield signals at the early LHC experiments. The slepton discov- 
ery plot in the mo - mi/2 plane (see [36] Fig. 13.29) shows that for rriQ = 120 GeV, our 
common choice for all three scenarios, the reach in mi/2 for 10 fb~^ is only about 150 GeV. 
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SIGNAL 


A 


B 


C 


a(pb) 


15.58 


5.74 


0.74 


2r + X 


10.9 


5.0 


1.2 


2e + X 


2.5 


2.8 


9.8 


2r + 06 + X 


11.5 


4.3 


1.3 


2r+(> lb)+X 


5.3 


2.9 


0.5 


2e + 0b + X 


2.8 


2.8 


14.8 


2e+(> 16) +X 


1.4 


1.8 


3.2 



Table 14: The S/VB{tt + QCD) ratio of the signals in Table 13 for C = lOfb' 



However, for (50 - 60) fb~^ sleptons in scenarios A and B should be detected. In the con- 
ventional r-annihilation scenario larger C will be required for this discovery. The direct 
detection of sleptons will provide additional insight into the DM relic density production. 

Several earlier analyses of DM relic density in niSUGRA or other models considered 
non-zero input values of Aq [37, 38] with different emphasis Arnowitt et al in [6] showed that 
the f-coannihilation corridor in the mo-''^i/2 plane is highly sensitive to Aq. For large Aq, 
f-coannihilation is effective even for rriQ as large as 1 TeV. In a more recent work [39] new 
bench mark points with non-zero Aq allowed by WMAP data were introduced. These points 
corresponds to new mass hierarchies and, consequently, new collider signatures. However, 
the correlation between leptonic signatures at LHC and different DM relic density producing 
mechanisms was not discussed before. 

In this paper we have considered the observed DM relic density constraints and direct 
constraints on sparticle masses from accelerators. We have not considered indirect constraints 
like the measured value of BR( B -^ 57). Many of the latter constraints arise from flavour 
violating processes. Strictly speaking these constraints are sensitive to the assumption that 
the quark and the squark mass matrices are aligned in the flavour space so that the same 
mixing matrix as the CKM matrix operate in the squark sector. This assumption of minimum 
flavour violation fails even if there are small off-diagonal elements of the squark mass matrix 
at the GUT scale. On the other hand such small elements does not affect processes like 
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neutralino annihilation and squark-gluino production and decay. For further discussions on 
this point we refer the reader to Djouadi et al in [21] and references there in. We, note in 
passing that model B and model C are allowed by the above constraint (see [18]). 

4 Summary and Conclusions 

We have examined the parameter space of the mSUGRA model with moderate values of 
the parameter tan (3. We focussed on zones of the rrio — 1711/2 plane corresponding to light 
sleptons and relatively light squark and gluinos compatible with the DM relic density data 
and constraints from direct sparticle searches. This part of the parameter space is interesting 
since viable signals from squark-gluino events are expected in the early stages of the LHC 
experiments. 

If one employs the often used but rather ad hoc assumption that the common trilinear 
coupling (Aq) vanishes at the GUT scale it is well known that there is only one such zone 
where r-coannihilation is the dominant mechanism for producing the relic density [21]. In 
this scenario the Xi and X2 present in squark gluino decay cascades exclusively decay via two 
body modes into appropriate neutrino (sneutrino) - slepton (lepton) pairs and lepton-slepton, 
neutrino-sneutrino pairs. Lepton flavour universality holds in these decays to a very good 
approximation. The sneutrinos decay into invisible channels with almost 100% BR. The X2 
also decays into invisible channels with large BRs. The suppression of hadronic decays of the 
lighter electroweak gauginos and the presence of additional carriers of missing energy lead 
to spectacular collider signals as has already been noted in the context of LEP/NLC [29, 30] 
and Tevatron [28]. In this paper we emphasize that this VLSP or ELSP scenario is realized 
in the popular r-coannihilation region of the mSUGRA model and discuss the signatures at 
the LHC in detail by introducing the benchmark scenario C (see Table 1). 

If the ad hoc assumption of vanishing Aq is given up additional WMAP allowed parameter 
spaces open up [18]. It is possible that the LSP pair annihilation and f- coannihilation 
both contribute significantly to DM relic density production although the former dominates. 
To illustrate the collider signals in this case the bench mark scenario A (see Table 1) is 
introduced. It is also possible that even with non-zero Aq the f coannihilation is the dominant 
mechanism but the corresponding squark-gluino masses are much smaller compared to C (see 
scenario B, Table 1). In scenario A and, to a lesser extent, in scenario B both x^ and X2 
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decay dominantly into final states involving r leptons (see the discussions in section 2 and 
Table 2) and lepton universality is violated. Moreover, these final states from squark-gluino 
production should be rich in 6-jets. 

In [18] the collider signals from squark-gluino events in three scenarios were compared 
and contrasted by exploiting the above characteristics. The calculations were done mostly 
at the parton level although in order to simulate the effects of r tagging some Pythia based 
analyses were also reported. 

In this paper we first employ the generic m-leptons + n-jets + ^t signatures to discrim- 
inate among the three scenarios ( see Tables 6 and 7), where only final states with stable 
leptons have been considered. We demonstrate that some qualitative idea about the DM relic 
density producing mechanisms may be obtained even without fiavour tagging. It is shown 
that the fraction of events with m = is much larger than that for m > 1 in scenarios A 
and B. The relative weight of the leptonic events, however, is significantly larger in scenario 
C (see Tables 5-7). 

Next we illustrate the r dominance of the final states in A and B by including r detection 
efficiency in our simulation (see Tables 11 -14). Here we have extended the analysis of [18] by 
considering 21 + X states where / stands for e and tagged r and X corresponds to hadronic 
states. In particular the observation of 2r -|- X events may provide a very convincing test of 
scenario A. The number of tagged 6-jets may further help to discriminate among different 
scenarios. Most of the crucial signatures discussed in this paper may be observed with C of 
10 fb~^. Some of them are observable with much smaller accumulated luminosity. 
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